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Abstract. We construct a perturbative expansion of the scalar sector in the Myers-Pospelov model, 
up to second order in the Lorentz violating parameter and taking into account its higher-order 
time derivative character. This expansion allows us to construct an hermitian positive-definite 
Hamiltonian which provides a correct basis for quantization. Demanding that the modified normal 
frequencies remain real requires the introduction of an upper bound in the magnitude |k| of the 
momentum, which is a manifestation of the effective character of the model. The free scalar 
propagator, including the corresponding modified dispersion relations, is also calculated to the 
given order, thus providing the starting point to consider radiative corrections when interactions 
are introduced. 
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1. INTRODUCTION 

The basic aim of this and previous work ^ is to define a consistent quantization of 
the Myers-Pospelov (MP) model considered as an effective field theory which provides 
perturbative corrections to standard QED. This model incorporates Lorentz invariance 
violating (LIV) corrections to QED, which are codified by dimension five operators 
which are C, T conserving but P violating. Our main concern with the construction, 
which we expect to play a major role in defining the appropriate way of making sense 
of the effective model, is that we recover the Lorentz invariant results of standard QED 
when such dimension five operators are turned off. In other words, we are looking for a 
smooth interpolation procedure between a Lorentz violating description and a Lorentz 
invariant one for QED. This point of view has been successfully carried on in Ref. [0], 
where LIV is codified by a dimensionless parameter modifying the integration measure 
in momentum space appearing in the calculation of one-loop processes. In this case LIV 
arises only from radiative corrections associated to standard particle Lagrangians; that 
is to say the zeroth order approximation is just the Standard Model of particles. The fact 
that LIV should be treated as a perturbation over the standard dynamics is supported 
by the very stringent limits which terrestrial experiments together with astrophysical 
observations set upon the parameters which label such violation. This requirement poses 
additional interesting challenges to this problem because we need to deal with higher- 
order time derivatives (HOTD) theories. The perturbative treatment of these, over the 
normal lower-order time derivative cases, is a non-trivial task, but fortunately it is well 
described in the literature [H, 3,(1, 0,0]. As an additional motivation to deal with quantum 
corrections in this model we mention the fine tuning problems recently reported dsj]. 



In this work we will concentrate in the scalar sector of the model, that can be used as 
the charged matter realization of scalar electrodynamics, which will be considered in a 
future work. Preliminary descriptions of the photonic and fermionic sector were already 
given in Ref. [lj]. We will set up the correct zeroth order quantization, identifying the new 
propagating normal modes described by their modified dispersion relations (MDR) up to 
second order in the LIV parameter. We will also present the corresponding propagator, 
with the final aim to calculate radiative corrections which are quadratical in the LIV 
parameters. 



2. PERTURBATIVE APPROACH TO HIGHER ORDER TIME 

DERIVATIVE THEORIES 

The general method for dealing with the canonical description of HOTD theories was 
given a long time ago in Ref. |9|]. In order to highlight some general features of these 
theories we briefly review their basic properties in a non-degenerate mechanical setting 
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7^ 0, where the fields depend only upon the time coordinate and 
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The generalization incorporating additional coordinates q(t) — > qi(t),I = l,...,N is 
direct and the introduction of space dependent fields goes along similar lines as in the 
standard transition from mechanics to field theory. 

If the highest time derivative in the Lagrangian is of order k , L = L(q(t) : ...,q( k \t)), 
the corresponding phase space will be of dimension 2k, been characterized 
by k coordinates: Qq — q(t) , Q\ = q^\t),...., Qk-i =q( k ~ 1 \t) together with k momenta 

m = ^+Z {-j t ) ^nj, i=o,...,(*-i). (2) 

In the sequel we avoid writing the explicit time dependence in the fields. The equation 
of motion for q will be of order 2k requiring the fixing of 2k initial conditions, which is 
consistent with the existence of 2k degrees of freedom in phase space. The Hamiltonian 
is 
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H(Qo,...,Q k -uPo,.^P k -i) = Y,P l Q l -L(Qo,---,Qk-i,q {k \Pk-i,Qo,...,Qk-i)) (3) 
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where the non-degeneracy assumption implies that we can solve q( k > as a function of 
Qoi •••) Qk-h Pk-l- The above expression is linear in the momenta P,, i = 0, ...,k — 2, thus 
making the Hamiltonian © unbounded from below, independently of the interaction 
terms included in the Lagrangian. 

Since we are interested in dealing with HOTD corrections in the action as perturba- 
tions upon standard theories, we must rely on a perturbation procedure which (i) retains 



the original number of degrees of freedom and (ii) produces a free Hamiltonian which 
is bounded from below as an adequate starting point for quantization. Such a method 
has been already developed in Refs. |@] and we present here a brief summary of it. 
To point out some of its basic features let us consider the simplest framework of a La- 
grangian depending upon accelerations L(q,q,q), where the HOTD contribution arising 
from jgq 2 is only present as a perturbation characterized by a small parameter g 

1 2 

L(q,q,q)^Lo(q,q) + ^gq ■ (4) 
The standard procedure of extremizing the action leads to 

8S=8^ 2 dtL(q,q,q) = [P Q 8q+P l 8q]f i + jT dtE (q,q,q,q (4) ) 8q(t), (5) 

where 

E(q,q,q,qW)=0, (6) 

is the fourth order equation of motion. From Eq. © the momenta can be directly read 
off as 

gq- (7) 
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dq 

From the simple form assumed for the HTOD term, which is appropriate for the non- 
degenerate situation, we can solve for q in terms of Pi and q^ in terms of Pq. Neverthe- 
less, notice that both substitutions carry the non-analytical factor l/g. This is precisely 
what makes non-trivial a perturbative expansion around g = 0. 

The Hamiltonian H and symplectic form Q. are defined according to the Ostrogradski 
method as 

H(q,q,P ,P l )=P q + Piq-L, £1 = dPo Adq + dP\ A dq. (8) 

The dangerous contributions to the Hamiltonian arise from the non-analytic term Pf /2g 
together with the unbounded piece Poq. 

Let us summarize now the general perturbation scheme for the non-degenerate case 
developed in Ref.[6], which applies to systems of the general form 

L = L Q (q,q)+gL l (q,q,...,q^). (9) 

In order to obtain the appropriate Hamiltonian to order g k , one starts by solving the 
equations of motion to order g^ -1 ) under the requirement of expressing all higher order 
derivatives , n > 2 only in terms of q and q, which provides the first step to define 
the effective variables of the problem. This will allow us to rewrite 

H = H(q,q), £1 = (0{q, q)dq A dq, (10) 

from where we can read the bracket {q, q}. Next we look for an invertible change of 
variables Q(q, q),P{q, q) such that 



Finally, the Hamiltonian H{Q,P) = H(q(Q,P) , q(Q,P)) together with the Poisson 
bracket {Q,P} = 1 define the physical effective system to the order considered. This 
Hamiltonian will be bounded from below provided that the initial one obtained from Lq 
is. The effective Lagrangian is given by L(<2, Q) =PQ — H and the corresponding Euler- 
Lagrange equations reproduce those of the original system to the order considered. Since 
this Lagrangian is first order in the time derivatives, the quantization is straightforward. 
A formal proof of self-consistency to all orders is provided in Ref.[6]. As clearly shown 
in one of the examples of such reference, the physical meaning of reducing the original 
phase space to that generated by q and q in the perturbative formulation is to suppress 
the excitation of high energy modes in a way consistent with the exact evolution and 
only to allow further excitations of the low energy modes already present in the zeroth 
order system. The method has been generalized to field theory in Ref. B7Q. 



3. THE COMPLEX SCALAR FIELD 

Before applying the method of Ref. [|6|] to our problem we need to identify the reduced 
phase space of this sector of the model which arises due to the presence of constraints, 
as it is shown in the next subsection. 



3.1. Reduced phase space dynamics 



We start from 



Xsadartf*,*) = "f (d 2 + m 2 ) (j) + ig(j)* (/1%) 3 0, 



(12) 



which describes the charged scalar field extension of the original MP Lagrangian. In the 
rest frame = (1,0) the equations of motion are 



§ = V 2 — m 2 + ig(j> 



(3) 



(13) 



together with its complex conjugate, which exhibit the HOTD character of the model. 
Inserting a plane wave in the above equation provides the MDR as solutions of a cubic 
equation. The condition for recovering the standard energy-momentum relation in the 
limit g — > requires that this cubic equation has three real solutions, one of which will 
be non-analytical in g, while the remaining two will be recovered in the perturbative 
approach to be used in the sequel. The above requirement demands 



a 2 F 2 ^ 4 
8 Ek < 27 s 



k| 2 + m 2 



(14) 



which sets an upper limit to |k| in a way analogous to the photonic and fermionic cases. 
For future comparison we write here the MDR obtained directly by solving the cubic 
equation to order g 2 



co±(k) =± 



E k ± X -gE 2 + 5 -g 2 El 



(15) 



The coordinate space variables of the model are 

(16) 

and our notation is 

^ ) = ^> * (1) =*, * (2) = #- ( 17 ) 

The presence of third-order time derivatives in the equation of motion (PTSl . instead of 
fourth order ones, precludes the existence of constraints in the Hamiltonian formalism. 
In order to deal with this construction we have found it convenient to perform one 
integration by parts in the dimension five operator appearing in Eq. (fT2~l) . This leads 
to 

j^ = <^ + 0*(V 2 -m 2 )0-^*f (18) 
The canonically conjugated momenta corresponding to the coordinates (fT6l) are 

n r = <p- i g $, n = 0* + i£0*, (19) 
= -igp, n^=o, (20) 

respectively. The Eqs. (l20l) are primary constraints and the extended Hamiltonian density 
is 

j?k = iv - - r (v 2 - m 2 )^> + -n r n^ + + p (m + i g <j)A . (21) 

The evolution of the primary constraints fixes the Lagrange multipliers C£,/3, so that 
no additional constraints arise and the original ones are second class. The number of 
coordinate degrees of freedom per space point is then i (2 x 4 — 2) = 3 which exceeds the 
usual two associated with the standard charged scalar field. We impose the constraints 

xi = n^, xi= n^+/^*, (22) 

strongly by introducing the corresponding Dirac brackets and verify that the resulting 
brackets for the reduced space variables 

(0,n ), (r,n r ), (fn^), (23) 

remain unchanged with respect to the original Poisson brackets, thus preserving the 
original symplectic structure. The reduced phase space Hamiltonian density is then 

Jtk = nrf + -Ilj (n r - <j>) - 0* (V 2 - m 2 )(j), (24) 

which exhibits a non-positive definite character due to terms linear in the coordinates 
and momenta; as well as the non-analytic behavior in l/g that prevents naive attempts of 
constructing perturbative schemes around g = 0. From Eq. (1241) we obtain the equations 



of motion 



n = (v 2 -m 2 )^*, n r = (v 2 -m 2 )<^, f^ = -(n -^n^, (25) 
i>* = l -n^, ^ = $ = - g (n r -<j>). (26) 

From them it is a direct matter to recover the Lagrangian equation (fT3l) and, consis- 
tently in an independent way, the corresponding complex conjugate expression. 



3.2. The perturbative expansion 

3.2.1. The symplectic form 

Regarding the symplectic structure of the theory let us indicate our conventions for 
the field theory case. Denoting generically Z A (f,x) the bosonic phase space fields (both 
coordinates and momenta), the Poisson brackets relations 

{z A (t,x),Z B (t,x')} = W AB (t;x,x > ), (27) 

are encoded in the corresponding symplectic two form €l(t) through the relation 

Q(t) = ^Jd 3 xd 3 x' W AB (t;x,x') dZ A {t,x) A dZ B {t,x'), (28) 

where Was(?;x,x') = — Wba(£;x',x), which can be considered as an antisymmetric ma- 
trix in its discrete as well as continuous indices, is such that 

j d 3 x' W A B(t;x,x , )W BC (t;x',x") = 5f S 3 (x-x"). (29) 

Our next step in the perturbative construction is to use the equations of motion to second 
order in g in order to rewrite the exact symplectic form 

Q= j d 3 x(dTl$ AdQ + dll r Ad^+dTl^ Ad<j>), (30) 

in terms only of the unperturbed variables (j), (j)* together with their first order time 
derivatives. 

The approximation to order g 2 leads to 

= (-E 2 ) [(1 + g 2 E 2 ) $ + ig<j)} + 0(g 3 ), (31) 

together with it complex conjugate, which we also substitute in the corresponding 
momentum expressions. Here we have introduced the notation 

£ 2 =(m 2 -V 2 ). (32) 



In this way the Hamiltonian density is 

JT = 0* (l - 2g 2 E 2 ) + ig^E 1 ^ - ig$*E 2 <p + 0*£ 2 0, 
together with the symplectic form 

a = J d 3 x(-2igd(j)*E 2 Ad(j)+d<j)* A(l-g 2 E z )d(j)) 
+ J d 3 x (d<p A (1 - g 2 E 2 ) d(j)* - igd<p* A d<j>) . 



(33) 



(34) 



From Eq. (1341) we read 
W AB (x,x') = 



+2igE 2 

-2igE 2 -(l-g 2 E 2 ) 

(\-g 2 E 2 ) 
(l-g 2 E 2 ) -ig 



(l-g 2 E 2 ) 


ig 




5 3 (x-x'). 
(35) 

The notation here is A,B = 1,2,3,4, with Z 1 = 0, Z 2 = 0*, Z 3 = and Z 4 = 0*. 
Inverting the matrix (T331) to second order in g we obtain 



<5 3 (x-x') 



(36) 

from where we read the non-zero brackets among the fields and their first order time 
derivatives. 








-ig 





l+3g 2 E 


W AB (x,x') = 


ig 





-(\+3g 2 E 2 ) 


l + 3g 2 E 2 





-2igE 2 




_-(l+3g 2 E 2 ) 





2igE 2 






{0(x,f),f (x^)} = -^5 3 (x-x'), 

{^(x,t),r^',t)} = {l+3g 2 E 2 )8 3 (x-x>), 

{ffrf)^^.*)} = {l+3g 2 E 2 )8 3 (x-x'), 

{<j)(x,t),<j)*(x',t)} = -2^ 2 S 3 (x-x'). 



(37) 
(38) 
(39) 
(40) 



3.2.2. New canonical variables 

The next step in the procedure is to introduce new canonical variables having the 
standard symplectic form in such a way that the final Hamiltonian density has the 
quadratic term in the momenta normalized to one. To this end we define, to order g 2 , 
the following new coordinates 0, 0* and momenta 11^, IT^* 

= [P(t> - f 2<A , = (R<j>*-i g £ 2 sr) , (4i) 



together with the corresponding complex conjugates: 0* = (0J , 11^* = (n^)*. The 
imposition of standard Poisson brackets among the tilde variables leads to the conditions 

P = Q, R = S, PR = 1- ^-E 2 . (42) 
The inverse functions are given by 



= ( 1+ i^)l^+^| n ^)- (43) 



<j) = (l + y^Viu^-UgE 2 A. (44) 
Substituting in the expression (1331 ) for the Hamiltonian density we obtain, to order g 2 , 

,w = n # n # . + 1 f>£ 2 % - n,#$) + r {e 2 + ^ (£ 2 ) 2 ) I (45) 

where the momenta are given by 

The Hamiltonian density can be conveniently written as 

je = W** +r (£>+t(£*f)i. (47) 

3.2.3. The modified dispersion relations 
Finally we construct the corresponding effective Lagrangian density 

leading to the equation of motion 

**to +lg pM>. + p +g 2 ^ Hg) = (49) 

Using a plane wave type solution one finds the MDR for positive frequencies to be 



g>i (k) = ( 1 + + - (gE k yj , ft)2 (k) = E k ^1 - -gE k + -(gE k ) z j . (50) 

The energy dependence in the corresponding plane wave are e~ ,COlt and e lC02 ' respec- 
tively. Let us emphasize two points at this level: (i) the MDR (|50l ) coincide with those 
derived from the exact equation in (TT5T) . (ii) the theory is charge conjugation invariant, 
which means that the field <j)(g) together with the charge conjugate field <j>c{g) — <!>*(— g) 
are solutions of the equation of motion (l49l) . as can be directly verified. 



3.3. The quantization 



It proceeds along the standard lines of the complex scalar field, the only difference 
being that the MDR are given by Eq.(l50l). In fact, we introduce the independent set of 
creation-annihilation operators 



a(k),a t (k / ) = 5 3 (k-k'), b(k),rf(ls!) = 5 3 (k-k') 
with standard vacuum 

fl (k)|0) = = fc(k)|0). 
In terms of these operators, we expand the complex field as 

d 3 k 



<j>(x) 



(51) 
(52) 

(53) 



together with its hermitian conjugate. The normalization factor £2(|k|) is determined by 
demanding the equal time commutation relation 



0(x),m(x') = /S 3 (x-x' 



and it is given by 



Cl{\k\)=E k [ l + ~-g 2 E 2 k 



5 

8< 



(54) 



(55) 



The remaining canonical commutation relations among the field operators can also be 
recovered, in virtue of the relations (|5TT) and (1551) . Starting from (|47l) . the normal ordered 
Hamiltonian turns out to be 



H 



0) 1 a t (k)a(k) + (0 2 ^ t (k)^(k) 



(56) 



This result can be considered as an additional consistency check of the choice (1551) 
to order g 2 . The interpretation of a' (k),a(k) Qr (k),£(k)) as creation and annihilation 
operators for positively (negatively) charged particles also follows directly. In fact, the 
normal ordered charge operator is given by 



Q 



a\k)a{k)-b\k)b{k) 



The Feynman propagator is 

iA F (x-y) = (0\T$(x)f(y))\0}, 
which can be separated into retarded and advanced pieces 

iA F (x-y) = e(x^-y )(0\^x)f(y)\0) + e(y -x )(0\f(y)^(x)\0). 



(57) 



(58) 



(59) 



The standard calculation yields 



^-y) = j ^q-^-EU^^}) TTe' m 



where the denominator reproduces exactly the momentum space version of the effective 
equation of motion ( l49i 
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